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ABSTRACT 


Predicting the effective life of materials for srace applicat;: r . - has become increasingly critical wiih the 
drive to reduce mission cost. Pr.'gram> have „ r.siderod nacy- - -muons to reduce launch costs auiudmg 
novel low mass materials and thin thermal blankets to reduce spacecraft mass. Determining the long-term 
survivability of these materials before launch is critical for mission success. This presentation will describe 
an analysis performed on the outer layer of the passive thermal control blanket ot the Hubble Space 
Telescope. This layer had degraded for unknown reasons during rhe mission, however ionizing radiation 
(IR) induced embrittlement was suspected. A methodology was developed which allowed direct 
comparison between the energy deposition of the natural environment and that of the laboratory generated 
env ironment. Commercial codes were u>ed t-. • predict the natural space IR environment, model energy 
deposition in the material from Null natural and laboratory [R N. urces. and design the most efficient test. 
Results were optimized for total and local energv deposition with an iterative spreadsheet. This method has 
been used successful!) for several laboratorv tests at the Marshall Space Flight Center. The stud) showed 
that the natural space IR environment, hv itself, did not cause the premature degradation observed in the 
thermal blanket. 
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1.0 INTRODUCTION 

The drive to reduce mission costs has resulted in innovative external materials for many applications in the 
space environment. Some of these materials have no relevant fight experience and are designed with very 
low tolerance for degradation. Thermal control blankets, which can be reduced in thickness tor a 
significant weight reduction, and mission specific optical films, which have little margin for degradation, 
are clear examples of potentially sensitive materials that are exposed to high ionizing radiation (IR) 
environments. The effect of the radiation environment, through mechanisms such as total ionizing dose 
(TID) and non-ionizing energy loss (NIEL), on sensitive materials like these impacts either the initial 
design or mission duration. 

The passive multi-layer insulation (MLI) thermal control blanket degraded significantly during the first 10 
years of the Hubble Space Telescope’s (HST) mission 1 . The MLI was inspected at 3.6 years and had 
developed minor cracks. After a total of 6.8 years on orbit, the blanket had developed severe cracks and 
was visibly degraded. Mechanical properties of samples returned to Goddard Space Flight Center (GSFC) 
were examined. This included tensile, bend, and micro-hardness tests. These tests indicated significant 
embrittlement, loss of ultimate tensile strength, and an increase in surface hardness; effects commonly 
observed in radiation induced embrittlement. There were other several natural environment sources that 
had impacted the blanket surface and each of those environments was tested independently to determine if 
there had been an individual source of the observed damage. Suspect environments that were tested 
independently included natural space IR. ultraviolet light, solar X-rays, atomic oxygen, and thermal 
cycling. 
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To determine the effects of the natural space IR environment on surface materials, the environment should 
be replicated as closely as possible in a laboratory experiment. The laboratory facility will be limited in the 
number or spectrum of energies available while the spectrum of the natural space environment is a 
continuum of energies starling from less than a keV and continuing through hundreds ot GeV. Since it will 
not be possible to replicate the natural space environment, some reference standard tor the eMect> of the 
natural environment must be determined and the laboratory environment will he generated to match me 
standard. The yardstick for comparison of the two environments was chosen to be total ionization energy 
deposition per unit mass (dose) as a function of depth in material. The units commonly used for dose are 
rads in material. 

While many forms of radiation induced material damage are not related to ionizing energy deposition, the 
ionizing damage is usually two or more orders of magnitude higher than the non-ionizing damage". This 
methodology could be used for non-ionizing dose from protons, but there are many additional code 
concerns with non-ionizing energy loss that were not considered here. 

There are several commercial and freeware codes that can model the natural space and laboratory 
environments and transport those environments through the spacecraft materials of intcre^. However, each 
of these codes has significant limitations, so that there is not a single code capable ot doing the entirety of 
this analysis. Each code is capable of performing some elements of the environmental analysis. By using 
the appropriate features of three separate codes and comparing their results using the reference standard, 
the complete analysis may be performed and verified 

Prediction of environment and initial dose as a function of depth was performed with the commercial code 
Space Radiation. This analysis provides spectra for the other code» and determines the relative significance 
of the proton and electron aspects of the ionizing radiation env ironment. 

The Integrated TIGER Series (ITS) is a Monte Carlo style code produced and updated by Oak Ridge 
National Lab. ITS was used to match the absorbed dose ot electrons which could be produced in the 
laboratory, to the absorbed dose of the natural space trapped electron spectrum on 10 mils (254 um) of FEP 
Teflon* for the orbit of the HST. The fact that the natural space environment will be an isotropic 
continuous spectrum, while the mono-energetic beams will be incident at 80° (near normal incidence) was 
considered along w ith the energy range limitations of the laboratory's electron beams. 

The Stopping and Range of Ions in Matter/Transport ot Radiation in Matter (SRIM 2000/TRIM) codes are 
ion transport Monte Carlo codes designed and updated by IBM Research. The SRIM 2000 codes were 
used to determine the effects of protons on the material tested. 

Space Radiation, using the internal code SHIELDOSE. was also used to verity the accuracy ot ITS and 
SRIM 2000. SHIELDOSE has poor resolution in the low energy range, but it has significant laboratory 
data on the energy deposition of high-energy electrons and ions . A comparison case which can be 
performed on Space Radiation and ITS is a thin plate of Aluminum with an isotropic electron input. To 
compare TRIM and Space Radiation, the range of hydrogen ions in a Hat plane of silicon was evaluated. 

Marshall Space Flight Center's Environmental Effects Group performed the tests. The laboratory used the 
test design resulting from the modeling to approximate the absorbed dose of the natural space IR 
environment for missions of one to ten years. The effects of those tests on the material samples were 
determined by measuring the optical properties. GSFC also performed tensile, bending, and micro-hardness 
tests on the samples. 

2.0 INITIAL PARAMETERS FOR ORBIT, ENVIRONMENT, AND MATERIAL. 

Damage to the MLI was assessed on a year-by-year basis for the 10-year reference mission for the HST 
with each year assumed to he equal to one tenth of the full mission. The approximate location of the HST 
over the time the sample was exposed was a circular orbit at 61 I kilometers and 28.5° inclination. 



The high-energy protons will also penetrate very deeply into a materia!. The average penetration depth for 
protons with energies greater than 10 MeV i.s greater than 10 mi A . From the plot of integral Hue nee in 
Figure F it can he estimated that less than 20# of the protons are of low enough energy to he stopped in 
the material 

The electron rluenee is significantly higher than the proton flue nee. For energies less than 500 keV, 

electrons have fluence more than 2 orders of magnitude greater than the protons. Because of these three 
factors, the electron fluence is the dominant source of the total dose in the first 200 mils of shielding. 

The margin of uncertainty for the fluence predicted by the AE-8 and AP-8 models is traditionally 

considered to be a factor of two*' 10 . Actual fluence of particles may vary by as much as three orders of 
magnitude during short periods. The radiation fluence is also non-isotropic so that the facing or direction 
of travel can influence the dose received and the large structures, such as the inner and outer Van Allen 
belts, also vary with time and specific solar conditions. These properties are not well modeled at present. 

It should he noted that the Jensen-Cain magnetic field model with epoch i960 was used throughout the 
period when trapped radiation data was collected for all AE/AP models’. These models should always be 
used with dates during the years in which data was taken. However, the magnetic fields have shifted in the 
last 40 years, so the predictions will be at different orbital locations than the spacecraft. There are versions 
of these codes being produced using current electron spectra, proton spectra, and magnetic field models, 
which will be good replacements. 

The Environmental Effects Group - Space Environment Effects Team at Marshall Space Flight Center 
performed the laboratory tests. In a vacuum chamber, proton and electron accelerators were both available 
providing mono-energetic beams of energies { .0 - 50.0 keV and 0.220 - 2.5 MeV. The beam lines, 
through which the radiation is emitted, are aligned at an inclination of 10 degrees from the normal to the 
test stand. The beams produced by the laboratory have a flux of approximately 6.25x10'* 
particles/cm Vsecond w ith even coverage of the test stand, which is square, three-inches on a side. 

The passive heat shield for the HST is a 10 mil (254 urn) layer of aluminized fluorinated ethylene 
propylene (FEP) Teflon \ The samples which were exposed in the laboratory were representative of the 
material used in orbit. The exact chemical composition and material preparation of this heat shield is not 
available, but the atomic composition and density is sufficient for radiation modeling. The material 
properties of Teflon* were given as 72.66# fluorine. 22.97# carbon, and 4.37 # oxygen with a density of 
2.15 grams/cm \ 

3.0 COMPARISON OF TRANSPORT MODELS FOR RADIATION IN MATTER 

3.1 Space Radiation with SHIELDOSE 

The commercial code Space Radiation uses the code SHIELDOSE to determine the energy deposition 
(dose) when transporting trapped proton and trapped electron spectra through a material. SHIELDOSE 
matches the ionizing radiation spectrum and target materials with pre-generated data 11 . There are two 
difficulties with SHIELDOSE in the Space Radiation format. The first is that the code uses pre-generated 
data, so that it doesn’t contain the information to generate new energy deposition profiles for new materials. 
Second, the code was not written for energies less than 100 keV. The code assumes particles of 100 keV 
are absorbed in the zone they are found, which creates an uncertainty of position on the scale of a 
micrometer. The stopping distances for higher energies are measured directly through experimentation, so 
that the error does not propagate throughout all of the results. 

Space Radiation uses any continuous energy spectrum of protons or electrons in units of 
fluence/crrr/second. This is propagated through the designated material with a result in dose (rads in 
material) as a function of depth in material integrated over mission duration. As Teflon* is not a material 
available in this package of Space Radiation, the initial comparisons of the four speetrums were made using 
silicon. 



The solar cycle also has a large effect on the radiation environment. During solar maximum, the electrons 
are at their peak values and the protons are at the lowest. During solar minimum, protons are at their 
highest and electrons are at their lowest fluence. 

While electron and proton cycles are opposite, the electrons are the usual driver for dose on external 
surfaces and protons are more significant for shielded materials. Since the material was recovered during 
solar minimum, it had been present during the worst of the solar cycle for external surfaces, solar 
maximum. The cycle lasts approximately 1 1 years between peaks. The solar output remains at maximum 
for only about a year and a half, but the entire 7 year period around the peak is modeled as the maximum. 
For the missions worst case radiation environment, only 3 of the 10 years are considered solar minimum. 

The trapped radiation environment was modeled using the Space Radiation commercial package, which 
includes AE-8 and AP-8. The environment was evaluated for the HST orbit using the AE-8 and AP-8 
codes for electron and proton fluence respectively. These codes each have a version for solar maximum 
and solar minimum^ 6 . The results for one year on orbit for AE-8 min. AE-8 max. AP-8 min, and AP-8 
max are presented in Figure 1. The galactic cosmic radiation and sokir event environments, such as those 
that occur during a coronal mass ejection or a solar proton event, were not included in the model because 
their effects are several orders of magnitude less than the effects of trapped ionizing radiation inside the 
Van Allen Belts 1 . 

The proton spectrum has a significantly greater energy range than the electron spectrum. The high energy 
protons have a much greater penetration depth than the electrons, but they do not have a high linear energy 
transfer rate until they are either captured or their total energy drops below 100 keV. So while each of the 
high-energy protons has more energy to deposit in a material, it will spread that energy over a much larger 
range than the electrons, effectively reducing the amount of energy deposited in a thin film. 
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Figure I Trapped Proton and Electron Spectra for Solar Maximum and Solar Minimum Conditions. 





In Figure 2, the dose per unit depth for each of the four spectrums is presented. It is clear that the ionization 
caused by electrons is two orders of magnitude greater than the ionization caused bv the protons in the first 
mil of the silicon and one order of magnitude throughout. More importantly, the electron dose has been 
shown to cause damage from as low as 10" rads in stressed fluoro-polymers. such as Teflon *’ In other 
recent experimental work on polymers, the material effects of low energy electrons, in high doses, have 
been clearly demonstrated 13 . Doses greater than Uf rads have been found to cause material damage in 
some materials, but doses above 10 rads cause material property degradation in most polymers 1 '. 

The proton dose, in comparison, is less significant by several orders of magnitude. Highly stressed 
electronic circuits and optical devices may also be effected by doses of 10 3 - 10 4 rads 14 , but no damage has 
been shown to polymers or other insulation materials in that range. 


Solar Minimum and Solar Maximum Annual Dose vs. Depth Curves 
Produce by SpaceRad for 61 1 km, 28.5 degrees Inclination 



Figure 2. Dose as a Function of Depth Curves in Silicon for the Natural Space Ionizing 
Radiation Environment for One Year of Orbit for the Hubble Space Telescope. 


3.2 Integrated TIGER Series 

The TIGER code was the Monte Carlo style code used to correlate the energy deposition results for trapped 
and laboratory electron environments. The electron spectrums presented in Figure 1 were used to represent 
the natural environment's isotropic spectrum. Single energy beams of the energies that could be produced 
in the laboratory at a 10° inclination from the surface normal were also modeled to represent the possible 
laboratory beams. The sample of Teflon* was divided into 100 regions in a thickness of 10 mils (0.0254 
cm) to measure energy deposition in each region. The low- energy cut off, in ITS. was reduced to 1 keV for 
greater accuracy when the step size was less than a few micrometers. This can be compared to 
SHIELDOSE, where the cut off cannot be adjusted and is set at 100 keV. For the regions very close to the 
surface, the low energy electrons will deposit the majority of the TID. 

The results of each TIGER run are presented in the units of MeV-crrr / gram-particle. To analyze these 
numbers, the units are divided into two parts as in equation 1 . 

MeV/g * cm' / particle ( 1 ) 


The first can be converted into rads. 






100 rads - 1 Joule / k 2 


( 2 ) 


1 MeV = 1.6021 XI O' ! ' ; Joules <3, 

so that the conversion factor can he found in equation 4. 

1 McV/g = 1 .602 1 X I O' rads (4) 

Equation 4 was used to convert the results to rads -cm 2 / particle. 

The second part of the equation is the inverse of the fluence. electrons/cm 2 , so that the results can be 
converted to rads for a given exposure. The natural environment has an integrated fluence for mission 
duration, while the laboratory environment can adjust the fluence produced so as to match to dose vs depth 
curves of the natural environment. The worst case of the electron spectrum in dose in Teflon* vs. depth in 
Teflon* is presented in Figure 3 with the factor of two uncertainty presented as error bars. 


Trapped Electron Dose vs. Depth in Teflon for 10 years at the 

Hubble Orbit 



Figure 3. The Dose vs. Depth Curve for 10 Years at the Orbit of the Hubble Space 
Telescope as Predicted by the Integrated Tiger Series 

The dose predicted in Figure 3 is significant for material properties. As stated previously, dose between 1 
and 10 Mrads causes degradation in some fiouro-polymers. The prediction does not, however, determine 
that that amount of dose will cause an effect in this sample, only that there is a reasonable cause for further 
examination of a specific material in the laboratory. 

3.3 The Stopping and Range of Ions in Matter/Transport of Radiation in Matter (SRIM 2000/TRIM) 

The SRIM 2000 model is a pure, Monte Carlo style, ion propagation code 16 . The code transports an 
individual mono-energetic ion in a layered and angled material using probability and path length to 
determine the ion’s energy loss to the material and the ion's total range. When the ion has stopped or passed 
through the material, the code repeats the process for the next ion. While an individual ion may produce 



any one of a number of energy loss interactions, the normalized sum of ten thousand ions usually converges 
to an accurate probability map. The number of ions used to generate the probability distribution can be 
adjusted to greater or lesser accuracy as required This process can be used to duplicate the laboratory 
environment, but cannot accurately reproduce the isotropic, spread spectrum of the natural environment. 

TRIM permits varying the ion species, energy and angle of incidence of the particle beam source and 
multiple layers with variation in thickness and material composition for each layer in the target. The 
output may be graphically represented as Monte Carlo style path tracks such as shown in Figure 4, or area 
integrated energy deposition in units of ionization energy per angstrom depth as seen in Figure 5. 



0 A - Target Depth - 254 um 


Figure 4. Monte Carlo Style Particle Traces Performed by TRIM Using a 
4 MeV Hydrogen Ion in 10 mils of Teflon* 

Figure 4 is a graphical representation of approximately 9.9xl0 4 ion paths produced by TRIM in 
determining the normalized ionization and energy deposition profiles presented in Figure 5. The primary 
information obtained from Figure 4 is that there is very little divergence between the paths. This leads to 
the conclusion that the number of ion paths used was sufficient to guarantee convergence. When there are 
a large number of high probability events, more ions may be needed. It is also clear that even 4 MeV ions 
did not penetrate the Tenon". As the laboratory environment was limited to 2.5 MeV, a much higher 
energy beam device would have been required to accurately test this sample. 

The results of the energy deposition profile are presented in units of eV/Angstrom/ion. To convert this to a 
unit of dose, the depth in material must be converted from angstroms, to units of mass, grams. This is done 
through a second step, converting the units of density from g/cm 3 to grams/cm 2 - Angstrom. Then, given a 
specific integral fluence, in units of ions/em*. the units can be converted to dose as follows. 

Normalized Energy Deposition (TRIM output) * Integral Fluence * 1/p = Energy Deposition/mass 

or 

eV/ Angstrom/ion * ions/ cm2 * (errr- Angstrom )/grams * l O' 6 MeV/eV = MeV / gram (5) 


The results of equation 5 can he converted, by equation 4. into rads in material. 

Figure 5 compares the probability of energy transfer for each method for each ion. Direct ionization from 
the impinging ion clearly deposits the majority of the energy. Direct recoil is very efficient in this energy 
range, but it is not a very common event until the ion has lost most of its energy. Recoil, the primary 
source of non-ionizing energy loss, generates approximately 2 orders of magnitude less energy than direct 
ionization. The recoiled ions themselves generate a small amount of ionization, roughly iQ7c of the total 
kinetic energy absorbed. Since Space Radiation has shown that the dose generated by protons should be 
approximately 1 kilorad per year, the non-ionizing dose should be approximately 10 rads/year. This is not 
a significant source of dose for materia! properties concerns. 



Figure 5. Normalized Energy Deposition Curves for Ionizing and Non-Ionizing Energy loss 

in Teflon e for 4 MeV Hydrogen Ions. 


3.4 Verification 

The codes were verified against each other by running comparison cases. A comparison that could be 
worked by both SHIELDOSE and TIGER was an isotropic exposure of a 10 mil thick aluminum plate to 
the solar minimum electron environment. The spectrum previously generated for solar minimum was used 
tor inputs to both TIGER and SHIELDOSE. The comparison spectrum is presented in Figure 6. The 
differences between the two results are obvious, but only a significant percentage of the result in the first 
mil (25.4 pm). The greater energy deposition, in the first few mils, from the TIGER series code represents 
its more detailed analysis of the low energy electrons. At deeper penetrations, lower energy electrons have 
lesser effects so that the codes produce nearly equal results. 



Comparison of Space Radiation and Tiger Series 



Figure 6. Comparison of Tiger Series Code against the Space Rad Code. 

A comparison that could be produced by both SHIELDOSE and TRIM was the penetration of hydrogen 
ions in Silicon. TRIM was run with single energy beams at normal inclination to the surface, SHIELDOSE 
was run with a stepped spectrum from which peaks emerged on the dose vs. depth graph. The results are 
presented in Table I . 


Energy 

SHIELDOSE 

TRIM 

100 keV 

< 1 micron 

0.8759 microns 

400 keV 

3.010 microns 

4.18 microns 

600 keV 

6.095 microns 

7.38 microns 

l MeV 

13.72 microns 

1 3.4 microns 


Table 1. Penetration depth of protons in Silicon as predicted by SRIM and SHIELDOSE. 

It is clear that SHIELDOSE and SRIM produce answers that differ on the order of a micron. However, as 
was seen in Figure 6, the values become less significant as the energy and distance into the material 
becomes larger. At 10 microns or IMeV ions, the results are less than 3 % different. These results are 
expected, SHIELDOSE was known to have this limitation, but Table 1 does show that the results from 
SHIELDOSE and SRIM are congruent at higher energies. 

4.0 RESULTS 

In Table 2 are the energies and total fluence used for matching the dose deposition of the beam and the 
natural space spectrum. For this result, the Chi square, £“» was 0.36, approximately a 1 -sigma divergence 
in a comparison data set. The energy deposition by these beams is approximately 1.4% less than the 
deposition by the natural space environment over 10 years. 

In Figure 9, the absorbed dose as a function of depth from the natural space IR spectra for the HST orbit is 
compared to the results from three of the electron beams. This is one example case of the combination of 
three beams. There are many variations in energy and tluence that could be used. By varying the initial 
choices in laboratory spectrum and different numbers of matching beams, the laboratory tests could be 
designed in several different ways. As there is a minimum time to turn the beams on and off of about one 
second, increased accuracy can be obtained by making the beam times as long as possible. The most 





5.0 CONCLUSIONS 


The methodology developed in this paper made it possible to accurately compare the natural space IR 
environment to the laboratory environment. The accuracy was very close, even considering the limitations 
of the test facility. The initial work with Space Radiation showed that electrons were the prime drivers of 
surface dose in low Earth orbit. ITS was able to model both the natural space and laboratory environments 
satisfactorily. While SRIM 2000 shows that the protons did not produce a level of energy deposition 
sufficient to cause damage to materials, the methodology, closely related to that used with ITS, is still 
valuable for other projects. All three codes are shown to correlate closely, which validates the methodology 
developed here. 

There are still areas where code improvements would be valuable. For example, there is not a model for 
isotropic proton irradiation for very thin materials. However, this function can be added to TRIM with a 
table containing isotropic entry angles and smooth variations in energy. The table would have to be 
generated by a Fortran code for each specific environment, but by using Monte Carlo processes there 
should be no reduction in final accuracy. There is also difficulty in measuring neutron production from 
certain reactions, such as Aluminums (p,2n) interaction with protons. Nuclear reactions are not modeled 
by TRIM, but neutron codes may be used for some predictions in the future. 

The laboratory tests showed no degradation of the sample for the ionizing radiation environment tested. 
This does not indicate that the electrons had no effect, but it does show that without contributory 
environmental or induced stresses, electrons were not responsible for the premature degradation observed 
on the HST MLI. The dose of electrons produced in this experiment was from 10 G to HT rads. The high 
surface dose caused by the low energy electrons may have had some effect, but there may not have been 
sufficient body damage to make a measurable optical or structural change. On the other hand, the electrons 
may not cause significant structural damage in Teflon® until greater than 10 s rads. There may be no 
damage done until much higher dose levels unless the material is stressed. 

There have been several tests with significantly higher doses, 10 9 rads, which showed electrons degrading a 
thin material of similar structure 12 . However, the test performed on those materials included simultaneous 
stresses with protons and UV radiation. The combined effects of the three sources of material stress may 
have had a synergistic effect that enhanced degradation. More testing will be required to determine the rate 
of material properties degradation of Teflon® from ionizing dose. 
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